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The present status of soft hadron production in high energy heavy-ion colhsions is 
•/^ I summarized. In spite of strong evidence for extensive dynamical evolution and collective 
expansion of the fireball before freeze-out I argue that its chemical composition is hardly 
changed by hadronic final state interactions. The measured hadron yields thus reflect the 
primordial conditions at hadronization. The observed production pattern is consistent 



\0 ■ with statistical hadronization at the Hagedorn temperature from a state of uncorrelated, 



color deconfined quarks and antiquarks, but requires non-trivial chemical evolution of the 
' fireball in a prehadronic (presumably QGP) stage before hadron formation. 
0^ ■ 

^ 1. Heavy-ion data and the nuclear phase diagram 

Relativistic heavy-ion collisions are studied with the goal of creating hot and dense 
y ! hadronic matter and to investigate the nuclear phase diagram at high temperatures 
pH I and densities, including the expected phase transition to a color deconfined quark-gluon 
plasma. But even if the energy deposited in the reaction zone is quickly randomized and 
the fireball constituents reach an approximate state of local thermal equilibrium, a simple 
connection between heavy-ion observables and the phase diagram is still not easy: the 
pressure generated by the thermalization process blows the fireball apart, causing a strong 
time dependence of its thermodynamic conditions which is difficult to unfold from the ex- 
perimental observations. There are therefore two fundamental issues to be solved before 
one can extract information on the nuclear phase diagram from heavy-ion experiments: 
(1) To what degree does the fireball approach local thermal equilibrium? (2) Which ob- 
servables are sensitive to which stage(s) of its dynamical evolution, and which is the most 
reliable procedure for extracting the corresponding thermodynamic information? 

Combining microscopic models for the dynamical fireball evolution with macroscopic 
thermal models for the analysis of heavy-ion data, significant progress has been recently 
made in answering both of these questions. Crucial for this achievement was the dra- 
matically improved quantity and quality of hadron production data from the analysis of 
collisions between very heavy nuclei (Au-|-Au, Pb-|-Pb) from SIS to SPS. 

Fig. 1 shows a compilation by Cleymans and Redlich [0 of hadronic freeze-out points in 
the nuclear phase diagram from various collision systems and beam energies. The upper 
set of points, parametrized by a constant average energy per particle (E) / {N) = l GeV [|I|, 
is obtained from measured hadron yields. They indicate the average thermodynamic con- 
ditions at chemical freeze-out when the hadron abundances stopped evolving. The lower 
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set of points, compared with lines of constant energy and particle density [0 , is obtained 
from analyses of hadron momentum spectra and/or two-particle momentum correlations. 
They indicate thermal freeze- out, i.e. the decoupling of the momentum distributions. The 
chemical and thermal freeze-out points at the SPS and AGS, respectively, are connected 
by isentropic expansion trajectories with S/B ^ 36-38 for the SPS and 12-14 for the AGS. 
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Figure 1. Compilation by Cleymans 
and Redlich [Q] of chemical and thermal 
freeze-out points. The legend refers to 
the symbols for the thermal freeze-out 
points. For the original references for 
all the data points see M. 
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My first goal is a critical discussion of how these freeze-out parameters were extracted 
from the data and how reliable Fig. 1 is. Following that is a more detailed study of the 
fireball properties at chemical freeze-out, taking into account additional information not 
contained in Fig. 1, and a discussion of a consistent dynamical picture which can explain 
Fig. 1. My main conclusion, based on a chain of arguments developed and sharpened 
over the last few years 0, is given in the abstract; similar conclusions were reached and 
recently publicized by R. Stock and are also found in E. Shuryak's talk H. 



2. Thermal freeze-out, "Hubble"-fiow, and the Little Bang 

Let me begin with a discussion of the thermal freeze-out points. Freeze-out marks the 
transition from a strongly coupled system, which evolves from one state of local thermal 
equilibrium to another, to a weakly coupled one of essentially free-streaming particles. If 
this transition happens quickly enough, the thermal momentum distributions (superim- 
posed by collective expansion flow) are frozen in, and the temperature and collective flow 
velocity at the transition "point" can be extracted from the measured momentum spectra. 
In high energy heavy-ion collisions the freeze-out process is triggered dynamically by the 
accelerating transverse expansion and the very rapid growth of the mean free paths as a 
result of the fast dilution of the matter |0. Idealizing the kinetic freeze-out process by a 
single point in the phase diagram is therefore not an entirely unreasonable procedure. 

As in the Big Bang, the observed momentum spectra mix the thermal information 
with the collective dynamics of the system. In the Big Bang, the observed microwave 
background radiation has a Bose-Einstein energy spectrum with an "effective tempera- 
ture" (inverse slope) which is redshifted by cosmological expansion down from the original 
freeze-out temperature of about 3000 K to an observed value of only 2.7 K. In the Little 
Bang, where we observe the thermal hadron radiation from the outside, the transverse mo- 
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mentum spectra are blueshifted by the collective transverse motion towards the observer. 
Simple approximate expressions which capture this effect are[| Tgiope ~ ^therm + ^rn{v±)'^ 



(which applies for p±<^m) and Tsiope ~ Tthermy i^+{'v±)) / {l—{v±)) (which is good for 
p_L3>m). For a given species (fixed m) the measured slope of the spectrum is thus am- 
biguous: temperature and flow cannot be separated. 

The ambiguity can be lifted in two ways: (i) One performs a simultaneous fit of the 
m_L-spectra of hadrons with different rest masses, thereby exploiting the mass depen- 
dence in the first of these two expressions 0. This makes the implicit assumption 
that thermal freeze-out happens simultaneously for all particle species. For Au+Au 
collisions at the AGS this works well and gives Ttherm ~ 93 MeV and {v±) ~ 0.5 c at 
midrapidity [§]. Or (ii) one concentrates on a single particle species and correlates (as 
described in detail by U. Wiedemann 0) their spectra with their two-particle Bose- 
Einstein correlations. The M^-dependence of the transverse HBT radius parameteif] 
R±{M±) ^ R/ ^JT+~^{v±y^^ then provides an orthogonal correlation between tem- 
perature and flow, allowing for their separation. For pions from Pb-|-Pb collisions at the 
SPS this leads again to kinetic freeze-out temperatures of 90-100 MeV and average tran- 
verse flow velocities of 0.5-0.55 c (perhaps even somewhat higher at midrapidity) |]9|,p!0[|. 

This fixes the position of the freeze-out point along the T-axis in Fig. 1, but what 
about /iB? Since chemical equilibrium is already broken earlier, at Tchem ~ 170 — 180 
MeV (see below), /^b is, strictly speaking, not well-defined. In order to still be able to 
associate with kinetic freeze-out a point in the phase-diagram one commonly adjusts fi^ 
in such a way that the deviations between the observed particle ratios and their chemical 
equilibrium values at Ttherm are minimized. This is acceptable if the deviations are small; 
in practice they can approach a factor 2 or so. This clearly causes irreducible systematic 
uncertainties in the baryon chemical potential at thermal freeze-out which are usually not 
evaluated and are not included in the horizontal error bars in Fig. 1 (where given). 

Let's nonetheless accept these thermal freeze-out parameters and now ask the question: 
How did the system get there? Does the implied picture of a rapidly expanding, locally 
thermalized fireball, the Little Bang, make sense? These questions can be studied by 



microscopic kinetic simulations of RQMD, URQMD and HSD |Tl| type; even if they 
do not include quark-gluon degrees of freedom during the very dense initial stages and 
thus may parametrize the initial hadron production incorrectly (see below), they can be 
used to explore the effects of scattering among the hadrons before kinetic freeze-out and 
the evolution of collective flow. A detailed study of thermalization by rescattering was 
recently performed by the URQMD group for Au+Au and Pb+Pb collisions from AGS 
to SPS energies [|I^] (Fig. 2). After an initial non-equihbrium stage lasting for about 8-10 
fm/c these systems reach a state of approximate local thermal equilibrium which expands 
and cools at roughly constant entropy for another 10 fm/c before decoupling. During the 
adiabatic expansion stage strong collective flow builds up. Thermalization is driven by 



^This is accurate for non-relativistic particles from a Gaussian source with a linear transverse velocity 
profile v^{r) = (t;j^)^r^, where {v±_) is the radial velocity at the rms radius r'^^^={x^+y'^). The analogous 
formula in lacks the factor ^ in the second term since it uses the radial velocity at r = rrms/V2- 
2^ ^ 0{\) accounts for different transverse density and flow profiles; i = \ for the case described in fn. [|. 
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intense elastic rescattering, dominated by resonances (e.g. vr+A^— i>A— >7r+A^); inelastic 
processes are much rarer and lead only to minor changes in the chemical composition of 
the fireball [1^,14]. As a result, significant deviations from chemical equilibrium occur 
which increase with time; most importantly, at thermal freeze-out one sees a large pion 
excess which can only partially be accounted for by the initial string fragmentation process 
Remarkably, these deviations from chemical equilibrium produce very little entropy 
The S/ B values extracted from the URQMD simulations agree with those from the 
thermal model analysis of the data (cf. Figs. 1 and 2). 
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Figure 2. Expansion trajecto- 
ries from URQMD simulations Wl 



Open and closed symbols denote the 
pre-equilibrium and hydrodynamic 
stages, respectively, of the collision 
in steps of 1 fm/c. The filled symbols 
lie on lines of constant entropy per 
baryon, 5/5=38, 20, 12 for 160, 40, 
10.7 A GeV, respectively The 
shaded region indicates the expected 
parameter range for the deconfining 
phase transition. 



Hadron momentum spectra and two-particle correlations thus provide strong evidence 
for the existence of the Little Bang: thermal hadron radiation with Ttherm ~ 90 — 100 
MeV and strong 3-dimensional ( "Hubble-like" ) expansion with transverse flow velocities 
^ 0.5 — 0.55 c (and even larger longitudinal ones These two observations play a 

similar role here as the discovery of Hubble expansion and the cosmic microwave radiation 
played for the Big Bang. But is there also a heavy-ion analogy to primordial Big Bang 
nucleosynthesis? In the following section I will argue that we have indeed evidence for 
"primordial hadrosynthesis" in the Little Bang. 



3. Thermal models for chemical freeze-out and "primordial hadrosynthesis" 

Chemical reactions, which exploit small inelastic fractions of the total cross section, 
are typically much slower than the (resonance dominated) elastic processes. One thus 
expects chemical freeze-out to occur before thermal freeze-out |T^ (Tchcm > ^thcrm) but on 
an expansion trajectory with roughly the same entropy per baryon S/ B. Fig. 1 suggests 
that this is indeed the case. This is analogous to the Big Bang where nucleosynthesis 
happened after about 3 minutes at Tchem ~ 100 keV whereas the microwave background 
decoupled much later, after about 300000 years at Ttherm ~ \ eV. The much smaller 
difference between the two decoupling temperatures in the Little Bang is mainly due to 
its much (about 18 orders of magnitude) faster expansion rate. 

Before discussing implications of the chemical freeze-out points in Fig. 1 I first explain 
how they were obtained. Can thermal models be used to analyze chemical freeze-out? 
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I discussed this question in some detail last year in Padova [0 and thus will be short 
here. The first difficulty arises from the collective expansion which strongly affects the 
shape of the m±- and ^/-spectra, in a way which depends on the particle mass. However, 
many experiments measure the particle yields only in small windows of m± and y. A 
chemical analysis of particle ratios from such experiments depends very strongly on model 
assumptions about the fireball dynamics W^ . Static ffieball fits yield chemical freeze-out 



parameters which are quite sensitive to the rapidity window covered by the data |]T7|. 
Flow effects drop out, however, from 47r-integrated particle ratios as long as freeze-out 
occurs at constant T and /i. An yields thus minimize the sensitivity to the collective 
ffieball dynamics and are preferable for thermal model analyses. 
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Of course, dynamic systems never freeze out at constant T and fi. While the steep T- 
dependence of the particle densities (which determine the local scattering rates and control 
freeze-out prohibits strong temperature variations across the freeze-out surface [|19 



incomplete baryon number stopping causes at higher energies significant longitudinal vari- 
ations of /iB [|16[- A global thermal ffi replaces the freeze-out distributions Tchem{x), fiB{x) 



by average values (Tchem), (/^b)- A recent study by Sollfrank in which he performed 
a global thermal ffi to particle yields from a hydrodynamic calculation, showed that after 
optimizing (Tchem) and (/xb) the thermal model predicted yields which differed by up to 
15% from the actual ones, although hydrodynamic simulations implement perfect (local) 
chemical equilibrium by construction. Thermal model fits can thus never be expected to be 
perfect] without detailed dynamical assumptions local variations of the thermal parame- 
ters in the real collision can never be fully absorbed by the model. Discrepancies between 
model and data at the 15-30% level are inside the systematic uncertainty band of the 
thermal model approach. While being impressed by how thermal models can reproduce 
the measured particle ratios at this level of accuracy over up to 3 orders of magnitude 
(see Fig. 3), I am thus deeply suspicious of "perfect" thermal model ffis. 



A second lesson to be learned from the exercise in ||20| is that the x^/d.o.f. resulting 



from such a fit is not very useful as an absolute measure for the quality of the fit: since 
discrepancies between the real yields and the predictions from the global thermal model 
cannot be avoided, x^/d-o.f. becomes larger and larger as the data become more and 
more accurate. While minimization can still be used to identify the optimal model 
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parameters within a given model, one should be very careful in using the absolute value 
of x^/d.o.f. to judge the relative quality of different model fits . 

Let me note that the ideal system for thermal model fits of particle ratios will be pro- 
vided by heavy-ion collisions at RHIC and LHC: hadron formation will happen at the 
confinement transition, and near midrapidity the baryon density is so low that is 
nearly independent of /xb- Replacing T{x) by (T) will then be an excellent approxima- 
tion. Due to longitudinal boost-invariance near midrapidity, knowledge of dN/ dy will be 
good enough for a reliable chemical analysis. Finally, transverse flow can only be stronger 
at RHIC and LHC than at the SPS, so freeze-out will happen even more quickly after 
hadronization, strengthening the primordial character of the observed particle ratios. 
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At the hadronization temperature 
Tc hadrons are created from the 
hydrodynamic phase with chem- 
ical equilibrium abundances and 
are then evolved kinetically by 
URQMD. 

Lower part: final hadron abun- 
dances at the end of the kinetic 
stage (circles) and if the calculation 
is stopped and all resonances are 
decayed directly at Tc (squares). 
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These expectations are borne out by a recent analysis by Bass and Dumitru [|l3l who 
combined a hydrodjTiamic description of the dense early stage with a URQMD simulation 
of the late hadronic stage. Fig. 4 (bottom) shows that indeed chemical freeze-out occurs 
quickly after hadronization: the yields at hadronization (squares) and after the last elastic 
scattering (circles) differ by less than 30%, in spite of many collisions in between |ll3[ . 

The fit of Pb+Pb data in Fig. 3 E^ yields a chemical freeze-out temperature Tchem ~ 
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170 MeV at full strangeness saturation (7s=l). In the same data are fit with Tchem 



144 ± 2 MeV and a strongly oversaturated strange phase-space (7s = 1.48 ± 0.08). The 
authors also allow for oversaturation of the light quarks and find 7q=1.72±0.08 which 
allows to absorb the large pion multiplicity at a low value of Tchem- This "chemical 



non-equilibrium" p2[ fit underpredicts Vt/'E by 40% and by 60%, a problem which 
disappears at Tchem=170 MeV (Fig. 3). More importantly, since 7^ = e^^l'^ , the freeze- 



out parameters in imply a very large pion chemical potential /i^r = 156 MeV > ttItt; 



this invalidates the assumption ||2^ that Bose statistics for pions can be accounted for by 
considering only the first correction to the Boltzmann term. Hence, while the authors of 
[p2| prefer their fit on the basis of a low x^/d.o.f., it has systematic uncertainties which 
far exceed the statistical errors given in ||22|| . 

Having established the location in the phase diagram where chemical freeze-out oc- 
curs, we should again ask: Row did the system get there? Since Tchem turns out to be 
very close to the predicted critical value for the hadronization phase transition, there is 
clearly no time between hadron formation and chemical freeze-out for kinetic equilibration 
of the hadron abundances by inelastic hadronic rescattering [Q. The observed hadronic 
chemical equilibrium at Tchem must therefore be pre-established: it reflects a statistical oc- 
cupation of the hadronic phase-space, following the principle of maximum entropy, by the 
hadronization process 0. Hadrons form from a prehadronic stage by filling each available 
phase-space cell with equal probability, subject only to the constraints of energy, baryon 
number and strangeness conservation (the latter includes a possible overall suppression 
of strangeness). Afterwards, the chemical composition decouples essentially immediately, 
without major modifications by hadronic rescattering. 

The parameter Tchem is thus not a hadronic temperature in the usual sense, i.e. not 
a result of hadronic kinetic equilibration. In the maximum entropy spirit it should be 
interpreted as a Lagrange multiplier which regulates the hadron abundances in accor- 
dance with the conservation laws and is directly related to the critical energy density at 
which hadronization can proceed. Tchem~170 MeV translates into ec~l GeV/fm^. This 



also explains naturally Becattini's observation p3|,|2[] of hadronic chemical equilibrium at 
the same value of Tchem in e~^e~ , pp and pp collisions at essentially all collision energies, 
although there hadronic final state interactions are completely absent. Whether the con- 
stituents of the prehadronic stage themselves thermalize before or during hadronization is 
a question which final hadron abundances cannot answer; as likely as their thermalization 
may appear, it is not necessary for an explanation of the observed phenomena. 

The concept of statistical hadron formation from a pre-existing state of color-deconfined, 
completely uncorrelated quarks and antiquarks is supported by a recent analysis^ of Bialas 
p5| which follows similar earlier arguments by Rafelski |2^ but formulates them more 
generally such that they do not require thermalization. Bialas points out that by consi- 
dering baryon/antibaryon (or generally particle/antiparticle) ratios, the unknown effects 
on hadronization from the internal hadron structure drop out and one can check directly 
whether the observed hadron abundances can be fully understood by just counting their 
conserved quantum numbers (carried by their valence quarks), or whether additional 
correlations among the quarks exist. He finds that in S-l-S and Pb-|-Pb collisions at 



^For an improved argument taking into account global flavor conservation during hadronization see 
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the SPS the former is true while hadron production yields in p+Pb collisions point to 
correlations among the quarks. 



4. Early memories: strangeness enhancement 

The one decisive feature which distinguishes heavy-ion from elementary particle colli- 
sions is the strangeness content in the hadronic final state: the global strangeness fraction 
of the produced quark- ant iquark pairs, Ag = 2(ss)/(M'U+(i(J)|produced, is about 2 times 
higher in nuclear collisions [p7| , p8| . This can not be reproduced by hadronic rescattering 



models p7| , p!^ and must thus be a feature of the prehadronic state. Here I would like to 
discuss in more detail the specific enhancement factors for K, A, A, S, S, Q, Q, and (j) 



reported recently and during this meeting [29,30 
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Fig. 5 shows that the relative strangeness enhancement between Pb+Pb and p+Be 
collisions is the stronger the more strange quarks the hadron contains. This is perfectly 
consistent with the above picture of statistical hadronization: disregarding other phase- 
space constraints, an Q, for example, which contains 3 strange quarks, is expected to be 
enhanced by a factor 2^=8 if strange quarks are enhanced by a factor 2. On the other 
hand, this pattern contradicts expectations from final state hadronic rescattering: since 
hadrons with more strange quarks are heavier and strangeness must be created in pairs, 
the production of stranger particles is suppressed by increasingly higher thresholds. 

An interesting observation from Fig. 5 is the apparent centrality independence of the 
specific strangeness enhancement factors: the enhancement appears to be already fully 
established in semiperipheral Pb+Pb collisions with about 100 participating nucleons. In 
fact, the global enhancement by a factor 2 was already seen in S+S collisions by NA35 



Since it must be a prehadronic feature, but the lifetime of the prehadronic stage is shorter 
for smaller collision systems, this points to a new fast strangeness production mechanism 
in the prehadronic stage. Exactly this was predicted for the quark-gluon plasma PT . 

At this meeting we saw new data on the centrality dependence of hadron yields [^]. Un- 
fortunately, different centrality measures and prescriptions to determine A^part have been 
used. This needs clarification before the pattern in Fig. 5 can be considered confirmed. 

Much recent effort went into trying to explain these observations within microscopic 
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simulations based on string breaking followed by hadronic rescattering. All such attempts 
failed. VENUS and RQMD give more strangeness enhancement for more central collisions 
p2 |, however not from hadronic rescattering, but mostly from the non-linear rise of the 



formation probability for quark matter droplets and color ropes. Strangeness enhance- 
ment is thus put in as an initial condition; unlike Fig. 5 it rises monotonically with A'^part- 
HIJING/i?i? uses baryon junction loops to enhance strange baryon production near 
midrapidity. Again this puts the enhancement into the initial conditions. The measured 
iVpart-dependence is not reproduced. The model also disagrees with the observed pattern 
fl/Q > S/S > A/A (which the statistical hadronization picture p5| explains nicely). 



Finally, the "improved dual parton model" which gets some fraction of the enhance 



ment in the initial state from "diquark-breaking collisions" and claims to obtain an even 
larger additional enhancement from hadronic final state interactions, suffers from a se- 
vere violation of detailed balance: it only includes inelastic channels (like ir+'B^Q+K) 
which increase multistrange baryons but neglects the (at least equally important [p5|j ) 
strangeness exchange processes (like n+Q^K+'E) which destroy them. Consequently it 
also fails to reproduce the apparent saturation of the enhancement factors seen in Fig. 5. 

The observed strangeness enhancement pattern thus cannot be generated by hadronic 
final state interactions, but must be put in at the beginning of the hadronic stage. No 
working model which does so in agreement with the data is known, except for conceptually 
simplest one, the statistical hadronization model. 

5. Summary 

The analysis of soft hadron production data at the SPS indicates that hadron for- 
mation proceeds by statistical hadronization from a prehadronic state of uncorrelated 
(color-deconfined) quarks. This leads to pre-established apparent chemical equilibrium 
among the formed hadrons at the confinement temperature T^; it is not caused by kinetic 
equilibration through hadronic rescattering. After hadronization the hadron abundances 
freeze out more or less immediately. The chemical freeze-out temperature thus coincides 
with the critical temperature, Tchem ~ Tc ~ 170-180 MeV, corresponding to a critical 
energy density ec ~ 1 GeV/fm^ as predicted by lattice QCD. 

The prehadronic state in {A>32) collisions contains about twice more strangeness 

than in e^e~ and pp collisions. This strangeness enhancement appears to be already fully 
established in nuclear collisions with 60 or more participant nucleons and can not be ge- 
nerated by hadronic final state interactions. This suggests a fast ss creation mechanism 
in the prehadronic stage, as predicted for a quark-gluon plasma. 

A clear hierarchy between chemical (Tchem ~ 170-180 MeV) and thermal (Ttherm ~ 90- 
100 MeV) freeze-out is observed in Pb+Pb collisions at the SPS; the gap is somewhat 
smaller (f» 130-140 MeV vs. ~ 90 MeV) at the AGS. In both cases thermal decoupling 
is accompanied by strong radial collective flow. The smaller inverse slopes of the Q m±- 
spectra suggest that a considerable fraction (but probably not all) of this flow is generated 
by strong elastic hadronic rescattering after hadronization 

I conclude that we have seen the Little Bang in the laboratory, and that most likely it 
is initiated by a quark-gluon plasma. 

Acknowledgements: I thank S. Bass, J. Cleymans and R. Lietava for providing me with 
Figs. 1, 4, and 5. Several fruitful discussions with R. Stock are gratefully acknowledged. 
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A comment by Zhangbu Xu led to the clarification presented in Footnote 1. 
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